Centrifugal separation and equilibration dynamics in an electron-antiproton plasma 
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Charges in cold, multiple-species, non-neutral plasmas separate radially by mass, forming 
centrifugally-separated states. Here, we report the first detailed measurements of such states 
in an electron-antiproton plasma, and the first observations of the separation dynamics in any 
centrifugally-separated system. While the observed equilibrium states are expected and in agree- 
ment with theory, the equilibration time is approximately constant over a wide range of parameters, 
a surprising and as yet unexplained result. Electron-antiproton plasmas play a crucial role in anti- 
hydrogen trapping experiments. 



PACS numbers: 52.27.Jt, 52.25.Fi, 36.10.-k 

Non-neutral (charged) plasmas held in a Penning- 
Malmberg trap rotate around the trap's magnetic axis. 
Such traps use a solenoidal field B to provide radial 
plasma confinement, and electrostatic fields to provide 
axial confinement. The rotation results from the E x B 
drift induced by the plasma's self electric field. If the 
plasma contains multiple species, the heavier will be 
pushed outwards by centrifugal forces. This effect, first 
predicted by O'Neil pQ, can lead to almost complete sep- 
aration of the species in sufficiently cold plasmas. Cen- 
trifugally separated states have been observed in several 
laser cooled Be + -ion [2 -4 and Be + -positron [5 systems. 

This paper presents the first images and detailed mea- 
surements of centrifugal separation in an electron (e - )- 
antiproton (p) plasma system. Figure [I] shows images of 
two such centrifugally-separated plasmas in the ALPHA 
antihydrogen trapping apparatus [6 j. (Recently, the 
ATRAP collaboration reported [7] indirect observations 
of centrifugally-separated states in a e~-p system.) We 
also report the first measurements of the separation dy- 
namics in any centrifugally-separated system. Specifi- 
cally, we report the timescale on which the p distribu- 



tion comes into equilibrium in response to changes in the 
e~ temperature and density. Depending on the type of 
change in the e~ parameters, the p equilibration time 
scale varies from milliseconds to seconds, but it is no- 
table that the time scale is approximately seventy mil- 
liseconds for a wide range of changes. Not all of these 
time scales can be explained by previously explored the- 
oretical mechanisms. 

The p 's will separate from the e~ 's when the E x B 
rotational velocity tcjr at the plasma edge starts to ex- 
ceed the thermal velocity of the p 's. In the limit, com- 
mon in this paper, that there are few p 's compared to 
e~'s, the p radial density will be given by np(r, z) ~ 
n e - (r, z) exp[(r/Lp) 2 ], where n e -(r,z) is the e~ radial 
density, and the p edge scale length Lp is 



" 2k h T ' 


1/2 


'8k h Te 2 B 2 ' 


1/2 






m^n 2 _e 2 

L P e J 





(1) 



where T is the plasma temperature, mp is the p mass, 
n e - = n e - (0, 0) is the central e~ density, and e is the 
unit charge. The rightmost expression holds only in the 
long plasma limit. 
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FIG. 1. (Color online) Images of centrifugally separated plas- 
mas trapped in a) a 1 T and b) a 3 T solenoidal field. In both 
cases the p plasma is the ring on the left, and the e~ plasma 
is the disk on the right. While the p and e~ plasmas are 
concentric in the trap, they do not image to the same center. 
The faint halos visible near the p plasmas are thought to be 
due to optical reflections. In these and subsequent figures, the 
color bar is linear in density (red highest) ; each image is indi- 
vidually adjusted to saturate the color bar. The 1 mm length 
bars give the true dimensions in the trap, not the dimensions 
on the MCP. 
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FIG. 2. (Color online) Centrifugal separation as a function 
of plasma size. In all cases the solid line is the e~ radial 
profile and the dotted line is the p radial profile. Apertures 
clip images D and E. 



Images like those shown in Fig. [T] are obtained de- 
structively [K,. One of the electrostatic well walls that 
confine the mixed plasmas is lowered, thereby allow- 
ing the plasmas to flow out of the uniform solenoidal 
region (IT or 3T) and onto a micro-channel plate 
(MCP) /phosphor screen/CCD imaging system, yielding 
an axially-integrated image of the plasma densities. The 
MCP /phosphor is located far into the solenoid's fring- 
ing field, in a field of only 0.024 T. Electrons follow the 
field lines well, but, because of their heavier mass, p 's 
are not as tightly bound to the field lines. We observe 
that initially overlapping e~ 's and p 's image to differ- 
ent locations; the locations depend on the originating 
magnetic field magnitude and the distance the particles 
travel to the MCP/Phosphor. When these two factors 
are held fixed, the centers of the e~ and p images are 
constant for all variations of particle number and radial 
profile, and are independent of the presence or absence 
of the other species. Further, the magnetic system is 
astigmatic for p 's, and the resulting images are always 
elliptical. The response of the imaging system to p 's is 
more than 100 x greater than it is for e _, s because the 
p annihilations greatly enhance the response [8], so the 
relative brightness of the e~ and p images can be deceiv- 
ing. For instance, in Fig. (TJd there are ~ 540 x as many 
e _, s than p 's. 

Figure [2] shows a scan of the equilibrium as a function 
of the radius r p of the e~ plasma. In all cases, the plasmas 
contain approximately 1.9 x 10 T e _, s and 35, 000 p 's, and 
the average temperature for the series is 120±30 K [9|fT0] . 
The radius of the plasma was controlled by varying the 
amount of radial compression [11 applied to the plasma 
before the equilibrium was established. The central e~ 
densities range from 1.9 x 10 9 cm -3 to 2.7 x 10 8 cm -3 . Un- 
der the constant total charge Q conditions of Fig. |2j_the 
relative scale ratio Lp/r p is proportional to Br p VT/Q. 
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FIG. 3. (Color online) Comparison of the measured radial 
profiles to computed radial profiles at bracketing tempera- 
tures. The experimental conditions are identical to profile 
Fig. [2]\. The solid lines are the nearly indistinguishable e~ 
radial profiles, and the dotted lines are the p radial profiles. 



Thus, we would expect that the e _, s and p 's would sep- 
arate less as the plasma radius increases, a trend clearly 
observed in Fig. [2] Our plasmas are close to ellipsoids 
with length to diameter ratios of 30; to more accurately 
estimate the degree of separation, we developed a nu- 
meric, two-dimensional (r-z) equilibrium code (N2dEC) 
which includes the necessary finite- length effects [2]-[4]. 
Figure [3] replots the smallest radius plasma from Fig. [2j 
and the results of N2dEC calculations at bracketing tem- 
peratures. As the measured temperatures are subject to 
systematic and shot-by-shot errors, it is not surprising 
that the best-fit temperature demanded by the code does 
not perfectly match the measured temperature. 

The line densities and (in later figures) moments 
used throughout this paper are computed from one- 
dimensional fits to the two dimensional images from the 
MCP/Phosphor. The fits also correct for the ellipticity of 
the p images; charge along all points on any given ellipse 
is mapped to the major radius of that ellipse. 

The e _, s in our plasmas cool via cyclotron radiation 
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FIG. 4. (Color online) Temperature and hollowness if as a 
function of time in a cooling plasma for a plasma similar to 
that in Fig. ^p. The predicted temperature is fit through the 
data at 0.25 s. 



with a calculated e-folding time of 0.427 s at 3T [12]; 
the e~ 's cool the p 's by collisions. If we actively heat 
a mixed e~-p plasma, we can watch the centrifugal sep- 
aration reemerge as the plasma cools. Figure [4] shows 
the measured and predicted temperature as a function of 
time in such an experiment. The measured temperatures 
follow the predicted temperatures until they level off, for 
unknown reasons, at about 130 K. Figure [4] also shows the 
measured hollowness H of the p radial profile as a func- 
tion of time, and H predicted by N2dEC calculations at 
the measured temperatures. This metric is defined by 
H = ( r 2 } - 5 /( r 4 }°- 25 ? where () denotes an average over 
the radial profile. H = 0.841 for a Gaussian profile, 0.915 
for a constant density ellipsoid, and 1 for a thin annulus. 
The value of H is robust to many forms of experimental 
noise, and does not depend on the overall plasma radius, 
but is subject to a small offset if the background in the 
images is imperfectly subtracted. Figure [4] shows that the 
p hollowness H follows the plasma temperature, suggest- 
ing that the p 's go through a series of quasi-equilibrium 
states. No delay in attaining equilibrium is visible. 

If we change the e~ conditions suddenly, the p 's must 
relax to the new equilibrium corresponding to these new 
conditions. Figure [5] shows a study of this process when 
a specified fraction of the electrons is suddenly removed. 
This is accomplished by lowering an electrostatic well 
sidewall for a time short (< 1 /xs) compared to the time 
over which the p 's can move, but long enough for the re- 
quired fraction of e~'s to escape. This process, called 
"e-kicking", substantially heats the remaining e _, s to 
several thousand Kelvin. Several cases were studied: 
leaving 10% of the e _, s, so that the density of e~'s is 
3.3 x 10 8 cm -3 ; leaving 0.5% of the e~'s, so that the den- 
sity is 5.3 x 10 7 cm -3 ; leaving 0.5% of the e~'s, but with 
3, 100 p 's rather than 35, 000 p 's as in all the other data; 
leaving 0.5% of the e~'s, but at 1 T rather than at 3 T as 
in all the other data; and leaving 0% of the e _, s. 

The e-kicking process leaves the p 's in their initial, 
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FIG. 5. (Color online) Hollowness as a function of time af- 
ter e-kicking leaves the stated percentages of electrons in the 
plasma. The state before e-kicking is similar to that shown 
in Fig. [TJd. Images a) correspond to the 0.5% series. The 
e~ plasma is the faint disk visible in the lower right of the im- 
ages (but remember that there are still more e~'s than p 's.) 
The apparent smallness of the e~ plasma is an artifact of the 
faintness of the image. Images b) correspond to the 0% series. 



ring-like distribution, but, because the e~ density is now 
low and the temperature now high, the final equilibrium 
profile is Gaussian-like. Considering the first four, non- 
0% cases only, the time to re-establish equilibrium (some- 
what arbitrarily set at H = 0.91) is approximately 60 ms 
(slightly longer for the low p case.) The images (Fig.^) 
for this class of relaxation show no sign of instability; the 
center appears to fill in gradually and uniformly. The 
last case, in which all the e _, s are removed, is very dif- 
ferent. The relaxation is approximately 10 x faster, and 
proceeds (Fig. [5)3) via a classic £ = 1 diocotron instability 
[T3] . Three-dimensional molecular dynamics simulations 
verify that the plasmas are indeed £ = 1 unstable un- 
der these conditions. In the other cases, leaving even as 
few as 95,000 e~'s (the 0.5% case) reduces the rotational 
shear to below the threshold necessary to drive the £ = 1 
instability. 

We have also studied the injection of a very small ra- 
dius p "slug" into a pre-existing, larger radius, cold e~ 
plasma. We then track two distinct measures: the expan- 
sion of the p slug out to the radius of the e _, s, and the 
subsequent hollowing of the p 's. Figure [6] shows two such 
series, corresponding to the injection of a standard load 
of p 's (35, 000) into a standard e~ plasma (Fig. [l]>like, 
but with an e~ density of 1.5 x 10 9 cm -3 ), and injec- 
tion of a reduced load of p 's (8, 800) into the same e~ 
plasma. The p 's begin to expand immediately, reaching 
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FIG. 6. (Color online) Radius and hollowness as a function 
of time after injection of 35,000 p 's (circles) and 8,800 p 's 
(triangles). The top set of images corresponds to the 35,000 
p 's series. 



a plateau after approximately 80 ms. Hollowing begins 
after about 20 ms, and plateaus shortly after the plasma 
stops expanding. Injecting the p slug heats the e~'s, par- 
ticularly for the standard p load where the increase in 
temperature is approximately 150 K. The late-time evo- 
lution (t > 200 ms) visible in the standard p load series 
is probably due to the e~ plasma radiating away the en- 
ergy brought in during the p injection. The reduced p 
load heats the e~ plasma less (approximately 50 K); con- 
sequently, it becomes hollow somewhat faster. 

In conclusion, the p 's in an initially hot e~ plasma 
appear to follow an evolving equilibrium as the e~ cool. 
The evolution time scale is set by the cooling timescale. 
When all the e~'s supporting a well-developed p ring 
equilibrium are removed (Fig. [5^), the ring collapses via 
an t = 1 diocotron instability on a time scale of a few 
milliseconds. In all other observed cases (Fig. [5|, in- 
cluding the injection of a p slug into a cold e~ plasma 
(Fig. |6|, equilibration takes on the order of 60-80 ms. 
These cases range in e~ density over a factor of nearly 
30, in e~ temperature over more than a factor of 10, in 
magnetic field over a factor of 3, and in p number over a 
factor of 11. The relevant collision frequencies, gyro radii, 
Debye lengths, etc. will vary over a proportionate range. 
The equilibration time scale for the higher density, lower 
temperature data shown in Fig. [6] is compatible with the 
predictions of particle diffusion/mobility theory [14 , but 
the data in Fig. [5] are not; the densities for these cases 
are too low, and the temperatures too high. Test particle 
calculations [15] predict equilibration times that are too 
slow by several orders of magnitude, even for the data 
shown in Fig. [6j A more exotic theory appears to be 
required, perhaps like that in Ref. [16 . 

The results presented here are relevant to several of the 
processes employed in antihydrogen (H) trapping exper- 



iments at CERN [6j[T7]. For example, e~'s are used to 
sympathetically cool the p 's, but are generally removed 
(e-kicked) before synthesis of H is initiated. Our mea- 
surements are the first detailed exploration of the radial 
and temporal dynamics of the p equilibration during this 
process. We have found it necessary to minimize the ra- 
dius of the p 's to enhance the probability of trapping [18] ; 
we do this by compressing the e~'s in a still mixed e~- 
p plasma, and relying on equilibration between the two 
species to compress the p 's [TTj. Compressing the e~'s 
too quickly leaves some p 's behind; the measurements 
here give a lower bound for the compression scale time. 
Finally, in some proposed mixing schemes [19 , the p 's 
and e _, s are left unseparated; the effects of centrifugal 
separation must be considered in these schemes. 
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